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Abstract 



We perform a quantitative analysis of the capability of K2K, MINOS, 
OPERA and a neutrino factory in a muon collider to discriminate the 
standard mass induced vacuum oscillation from the pure decoherence 
solution to the atmospheric neutrino problem and thereby contribute 
to unravel the dynamics that governs the observed Vn disappearance. 
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I. INTRODUCTION 



Two years ago, the Super-Kamiokande (SuperK) atmospheric neutrino data as- 
tonished the world giving the first compelling evidence in favor of — > v T oscil- 
lation This incredible result has since been confirmed by other atmospheric 
neutrino experiments PH, as well as by the preliminary K2K disappearance 
experiment result H, making unquestionable the fact that neutrinos suffer flavor 
conversion. 

Naively, one may think that proves neutrinos have non-zero mass and that the 
next challenge for experimentalists is simply to determine the neutrino mass squared 
differences and the texture of the neutrino mixing matrix. Indeed, if the dynamics 
of neutrino flavor change is mass induced in the standard way H , this is obviously 
the next logical step. Unfortunately, this is not an established fact. 

Although the atmospheric neutrino data collected up to now allow one to definite 
exclude some energy dependencies for the — > u T conversion probability ||, some 
interesting possibilities, such as neutrino decay [0] and pure quantum decoherence 
are capable of explaining the data comparably well to the standard mass induced 
oscillation mechanism. This in spite of the fact that the dynamics behind neutrino 
decay and pure decoherence give rise to a survival probability monotonically 
decreasing with the neutrino energy while the mass induced mechanism leads to 
a harmonic probability of oscillation. We therefore believe that another important 
experimental task should be to unravel the nature of the flavor changing mechanism. 
We have to point out that the neutrino decay scenario is also mass induced, but 
throughout this paper when we allude to the mass induced mechanism, we will be 
referring to the standard neutrino flavor oscillation scenario. 

Many, if not all, of the proposed future neutrino long baseline experiments were 
designed to measure — > v T oscillations in order to pin down the oscillation pa- 
rameters having in mind the standard mass induced oscillation mechanism. It is 
important to verify their real capabilities to discriminate among different flavor 
changing dynamics. 

The purpose of this paper is to investigate to what extent K2K and the next 
generation neutrino oscillation experiments will be able to discriminate the mass 
induced — > v T oscillation solution || to the atmospheric neutrino problem ( ANP) , 
from the recently proposed pure decoherence one |8[]. We do not investigate in this 
paper the decay mechanism, since it implies in the existence of sterile neutrinos, 



giving rise to a richer phenomenology, this study will be reported elsewhere \W . 

The outline of the paper is as follows. In Sec. |J, we briefly review the mass 
induced and the pure decoherence mechanisms of flavor conversion. In Sec. [TI], we 
define the statistical significance tests we will use in order to quantify the separation 
between the two ANP solutions. In Sec. [TV], we discuss the power of discrimination 
of the mass induced oscillation solution to the ANP from the dissipative one at 
K2K fill, MINOS pi, OPERA pf and a possible neutrino factory in a muon 
collider [14]. Finally, in Sec. [V], we present our conclusions. 
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II. REVIEW OF THE FORMALISM 



The time evolution of neutrinos created at a given flavor by weak interactions, 
as of any quantum state, can be described using the density matrix formalism by 
the Liouville equation [PT5| . If we add an extra term L[p^\ to the Liouville equation, 
quantum states can develop dissipation and irreversibility fl5l , |16l 
Liouville equation for p M (t) can then be written as [|15[] 



The generalized 



dt 



-i[H, PlM (t)] + L[ Ptl (t)}, 



where the effective hamiltonian H is, in vacuum, given by 



H 



A 
-A 



(2) 



where A = (m 2 — mf)/AE l ,, we have already considered ultra-relativistic neutrinos 
of energy E u and the irrelevant global phase has been subtracted out. We assume 
here oscillation only between and v T in a two generation scheme. 

The most general parametrization for L[p] contains six real parameters which 
are not independent if one assumes the complete positivity condition |15] . In one of 
the simplest situation, which in fact physically arises when the weak coupling limit 
condition is satisfied, only one of the parameter, 7, has to be considered. In this 
limit, Eq. ([!]) can be solved to calculate the [|1^ 



P{vp ~^v T ) = - sin 2 29 [1 - e^ L cos(2AL)], 



(3) 



the probability of finding the neutrino produced in the flavor state in the flavor 
state v T after traveling a distance L under the influence of quantum dissipation 
driven by the parameter 7. 

When 7 = 0, we get the usual mass induced oscillation (MIO) probability in two 
generation. The survival probability, in this case, is the standard one 



P( Vlx -». Vr ) = - sin 2 29 [1 - cos(2AL)]. 



(4) 



On the other hand, if neutrinos are massless or degenerate (Am 2 = 0) and 
weak interaction eigenstates are equal to mass eigenstates, even though standard 
oscillations cannot occur, flavor conversion can still take place through the pure de- 
coherence mechanism (PDM) |T7| . Explicitly, the neutrino flavor change probability, 
in the simplest case where a single decoherence parameter is considered, becomes 



1 



[1-e 



-2 7 L! 



(5) 



We will assume here for PDM that 7 = •y (E u /GeV)~ 1 , where 70 is a constant 
given in GeV. This Ansatz may be motivated by the assumption that the exponent 
in Eq. (EI) behaves like a scalar under Lorentz transformations |§. 
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Note that here we will be studying both MIO and PDM in a two generation 
framework since this is enough to explain well the atmospheric neutrino data. 
Notwithstanding, one may wonder about the contribution from electron neutrinos. 
If PDM really takes place in nature, it should in fact involve all three neutrino fla- 
vors. We can have some insight on what can happen in this case. From the results of 



Refs. pjnj on the limits on PDM from SN 1987A data, and with the Ansatz above, 



we expect the decoherence parameter that accompany the — > v e conversion to be 
smaller than 10~ 39 GeV. This means that for the energies and distances of the long 
baseline experiments we will study here, this effect should be completely negligible. 
Nevertheless, a complete study of the PDM with three neutrino generations should 
in fact be performed to confirm this assumption. 



III. STATISTICAL SIGNIFICANCE TEST 

In order to define the capability of an experiment to discriminate mass induced 
Vp — > v T oscillation from pure decoherence, we define the number of standard devi- 
ations of separation between MIO and PDM as n a = where 



X%o, sin 2 29, Am 2 ) = 2[iV PDM ( 7o ) - jV MIO (sin 2 29, Am 




(6) 



is the confidence level according to the procedure proposed by the Particle Data 
Book jnj. Here, -/V PDM ( 7 o) is the total number of events theoretically expected if 



PDM is the solution to the ANP and iV MIO (sin 2 29, Am 2 ) is the total number of 
events that can be observed by the experiment as a function of the two parameters 
involved in the MIO mechanism. 

For each experiment we have studied, we have computed two different types of 
contour level curves. 

First, we fix jV PDM = iV PDM (7Q est ), at the number corresponding to 70 = 7g est = 
0.6 x 10~ 21 GeV, the best-fit point of the PDM solution to the ANP g] and vary 
the MIO parameters in the interval 1 x 10~ 3 eV 2 < Am 2 < 2 x 1CT 2 eV 2 and 
0.8 < sin 2 29 < 1 consistent with the atmospheric neutrino data. In this way, we 
obtain curves of fixed n a in the plane sin 2 29 x Am 2 . 

Second, we fix sin 2 29 = 1 and vary the PDM parameter 7o in the interval 
0.25 x 10" 21 GeV < 70 < 1.1 x 10~ 21 GeV, as well as the MIO parameter Am 2 in 
the interval 1 x 10 -3 eV 2 < Am 2 < 2 x 10~ 2 eV 2 . The upper limit of the range 



in 70 is the one allowed by the CHORUS/NOMAD data |5,17|, the lower limit was 
estimated by using 70 ~ 2.54 x 10~ 19 x ( Am 2 /eV 2 ) GeV, with Am 2 = 1.0 x 10 -3 eV 2 . 
Like this, we can get curves of fixed n a in the plane 70 x Am 2 . This allows us to 
extend our conclusions beyond the best-fit value of the PDM solution to the ANP. 
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IV. PDM VERSUS MIO 



We have investigated the capability of K2K [ TT], 20|, 2 1 1 and the next generation 



neutrino oscillation experiments MINOS [O, OPERA [IT3|,P3] and a neutrino factory 



in a muon collider ||14|| , to discriminate the PDM solution to the ANP with 7 ~ 
0.6 x 10~ 21 GeV Q, using the ansatz 7 = 'j (E u /GeV)~ l , given in Ref. §|, from the 
traditional one due to — > v T MIO in vacuum with Am 2 ~ (1.1 — 7.8) x 10~ 3 eV 2 
and sin 2 29 > 0.84 ||. We recall that the best-fit point for the MIO solution to the 
ANP is at (sin 2 29, Am 2 ) = (1.0,3.0 x 10~ 3 eV 2 ) §. 

We would like to point out that the decoherence solution to the ANP is open 
to two different readings: either it can be viewed as an effect of pure decoherence 
or as a combination of quantum decoherence plus vacuum oscillation driven by 
Am 2 < 10~ 6 eV 2 with sin 2 29 ~ 1. In the first case, there is a single free parameter 
(70) and the flavor change probability is given by Eq. ((5|), in the second, there are 
two free parameters (70 and sin 2 29) since the probability will be given by Eq. (||]) 
with cos(2AL) -> 1. 

We now present and discuss the results of our study. 

A. K2K 

In Ref. ||, K2K was cited as a possible experiment to test the novel decoher- 
ence solution to the ANP. This possibility would be very appealing for K2K is an 
experiment which is currently taking data. 

In order to verify this, we have calculated the expected number of events in K2K 
for the goal of the experiment, i.e. 10 20 protons on target (POT) Q, for three 
hypotheses: no flavor conversion, mass induced v ^ — > v T oscillation with parameters 
consistent with SuperK atmospheric neutrino results || and the pure decoherence 
solution to the ANP @. 

K2K is a — > disappearance experiment, where the muon neutrinos that have 
an average energy of 1.4 GeV are produced by the KEK accelerator, first measured 
after traveling 300 m by the Near Detector, which is a 1 kton water Cherenkov 
detector, and finally measured by the Far Detector (FD), the SuperK 22.5 kton 
water Cherenkov detector localized at 250 km from the target. This experiment, 
which has started taking data last year and has currently accumulated 2.29 x 10 19 
POT, seems to be confirming — > disappearance as expected by the atmospheric 
neutrino results [|]. 

The expected number of events in the FD can be computed as follows 

N FD = RN^°, (7) 



1 We remark that in our notation 270 corresponds to 70 of Ref. 
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where R is the ratio between the number of observed over the number of expected 
events in the Near Detector, we have used the fact that R ~ 0.84 from Ref. ||20|| , 
and N F ^° is the theoretical expectation that can be calculated as 

= vn FD J ^ FD (E)a(E)P(y li -> uJdE, (8) 

where E is the neutrino energy, &fd(E) is the z/ M flux distribution at the Far De- 



tector, cr(E) is the total neutrino interaction cross section taken from Ref. [23| and 
ripD is the number of active targets in the FD. Also, we have introduced a normal- 
ization factor r) which was fixed to 0.65 in order to get the same expected number 
of events as K2K for null oscillation. The shape of <$>ft>(E), was taken from Ref. 



11|], but the total flux has been renormalized to account for the number quoted in 
Table 1 of Ref. [20|| . The survival probability P{p IJj — > v^) = 1 — Piy^ — > u T ) in two 
generations, with Pivp, — > v T ) either equal to zero (for no — > v T conversion), to 
the usual two generation MIO probability (Eq. (f|)) or the PDM flavor conversion 
probability (Eq. (|)). 

We first perform the calculation of the total number of expected events, iVpo, 
in the absence of any flavor change and for oscillation with sin 2 29 = 1 and Am 2 = 
3x 10~ 3 eV 2 , 5x HT 3 eV 2 and 7x HT 3 eV 2 for 2.29 xlO 19 POT. These results, which 
agree quite well with the K2K estimations presented in Ref. ||], are summarized in 
Table |. Thus we are confident that our numbers are reasonable and we can proceed 
to estimate the total number of events expected in the FD when K2K reaches 10 20 
POT for vacuum oscillation and decoherence. These numbers are also reported in 
Table |. From this, we can see that the mass induced oscillation and the decoherence 
effect at their best-fit values imply, for the goal of the K2K experiment, values for the 
total number of events which are statistically compatible, hence the two solutions 
will be indistinct at K2K. 

One may wonder about the energy distribution of the K2K events, 
which, in principle, could be used to discriminate the solutions. We 
show, in Fig. [I], the Monte Carlo simulated reconstructed neutrino spec- 



trum for one ring /i-like events at SuperK taken from Ref. |^TJ for Am 2 



0.0015 eV 2 , 0.0028 eV 2 , 0.005 eV 2 and 0.01 eV 2 and sin 2 29 = 1, as well as our esti- 
mation of the distortion expected for the best-fit point of the PDM solution to the 
ANP. The latter was done simply by multiplying the bin content for no oscillation 
by the average PDM survival probability in the bin. 

From Fig. [TJ, we see that if one takes into account only the statistic error, com- 
pletely disregarding the systematic one, the curves are already virtually indistin- 
guishable so that one unfortunately can hardly hope to discriminate between these 
two solutions with the K2K data. It is important to remark that the actual spec- 
trum is currently under study, so what we present here can be viewed as a tendency 
which indicates that, even if one compares only the shape for the best-fit values of 
the two solutions, this type of discrimination will be very difficult in K2K. 

In Figs. ^| and |3|, we show the result of the statistical significance tests for K2K, 



as proposed in Sec. III. These curves were calculated with the total number of 
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events N FD and, although this is not directly related to what is plotted in Fig. [I, 
it takes us basically to the same conclusion, i.e. for data compatible with 2.2 x 
10 -3 eV 2 < Am 2 < 4.0 x 10~ 3 eV 2 the maximal separation between the two solutions 
cannot exceed ~ 2 a if 70 = 7o est . We see in Fig. || that for data consistent with 
2.2 x 1CT 3 eV 2 < Am 2 < 3.5 x 10~ 3 eV 2 the separation is always less than 3 a, V70. 
A separation of 5 or more a can only be achieved if the data is compatible with 
Am 2 > 7 x 1CT 3 eV 2 . We also see that the point which correspond to the best-fit of 
the MIO and the PDM solutions is inside the n a < 1 region. Therefore, we conclude 
that it will be rather difficult to disentangle the two ANP solutions before the arrival 
of the next generation neutrino experiments. 



B. MINOS 



The MINOS experiment [[L2| is part of the Fermilab NuMI Project. The neutrinos 
which constitute the MINOS beam will be the result of the decay of pions and kaons 
that will be produced by the 120 GeV proton high intensity beam extracted from 
the Fermilab Main Injector. There will be two MINOS detectors, one located at 
Fermilab (the near detector) and another located in the Soudan mine in Minnesota, 
about 732 km away (the far detector of 5.4 kton). 

According to Ref. P3| , MINOS will be able to measure independently the rates 
and the energy spectra for muonless (0/i) and single muon (1/z) events, which are 
related to the neutral current (nc) and charged current (cc) reactions. Three different 
neutrino energy regions are possible: low [E v ~ 3 GeV), medium [E v ~ 7 GeV) 
and high [E v ~ 15 GeV). For MINOS to operate as a v T appearance experiment 
either the high or medium energy beam is required, since one must be above the r 
threshold of 3.1 GeV. 

We have studied here two different observables that can be measured in MINOS 
: the 1/x-event energy spectrum and the 0/i/l/z event ratio. The three different beam 
possibilities were investigated. 

The expected number of 1/x-events in MINOS, dNi^, can be calculated using 

= lt iK) {[1 " ~> Ut)] + P(l/ " Vt) Vm Br(T "> M ' (9) 

where rj(E u ) = a UT ^ cc (E u ) / a u ^ cc (E u ) , is the ratio of the cc-cross section for v T over 
the cc-cross section for z/ M , Br(r — > /i) is the branching ratio of the tau leptonic decay 
to muon (18%), dN cc /dE u is the energy spectrum for cc-events in the MINOS 



far detector in the case of no flavor change [24J] and Piy^ — > v T ) is the probability 
of — > v T conversion. 

There are two different ways a muon can be produced; either by the surviving 
Ufj, which interact with the detector (first term of Eq. (Bp) or by the contribution 
from taus generated by v T interactions in the detector, after — > v T conversion, 
followed by r — > u^fi decay (second term of Eq. @). In both cases, the events 
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must trigger the detector and be identified as muons to count as 1 //-events. Here, 
we have not considered the possible contamination from neutral current events and 
the trigger and identification efficiencies were supposed to be 100%. 

In Fig. ^, we present the 1/t-event energy distribution expected at MINOS for 
no flavor change and the best-fit parameters of the MIO and of the PDM solutions 
to the ANP for 10 kton year (~ 2 years of running). We can observe that if the 
final choice of beam is the low energy configuration, MINOS will be in the same 
footing of K2K, meaning that discrimination between solutions will be extremely 
disfavored due to low statistics. In the event of the choice fall on the medium or the 
high energy beam, discrimination will become more likely due to higher statistics. 

We show in Fig. [j| the result of the statistical significance test in the Am 2 x sin 2 29 
plane, for the three beam setups. We see that for Am 2 = 3 x 10~ 3 eV 2 , PDM and 
MIO are separated by less than 2 a (low), more than 3 a (medium) and more 
than 5 a (high), Vsin 2 2# in the range compatible with SuperK atmospheric data, 
confirming our previous conclusions on Fig. |j. 

We can further point out that, considering the 99 % C.L. sin 2 29 x Am 2 region 
allowed by the SuperK data ||, the two solutions will be indistinguishable within 3 
a for 1.7 x 10~ 3 eV 2 < Am 2 < 5.3 x 10~ 3 eV 2 in the low energy configuration, for 
3.0 x 10~ 3 eV 2 < Am 2 < 4.3 x 10~ 3 eV 2 in the medium energy configuration and 
for 4.2 x 10 -3 eV 2 < Am 2 < 6.2 x 10~ 3 eV 2 in the high energy configuration. 

In Fig. H|, we show the statistical significance test in the plane 70 x Am 2 for the 
three energy possibilities. The point which correspond to the best-fit of the MIO 
and the PDM solutions is inside the n a < 1 region for low, the n a < 2 region for 
medium and only for the high energy case it falls into the n a > 5 region. Here it is 
again demonstrated that the medium and high energy setups are preferable to the 
low one if the issue is to discriminate the two mechanisms. 

We have also computed the ratio Ro^/i^ that should be expected in MINOS as a 
function of the free parameters of the flavor changing hypothesis. This ratio has the 
advantage that it does not require the understanding of the relative fluxes at the 
near and far detectors, it is also quite sensitive to neutrino flavor conversion since 
when they occur not only 1/t-events are depleted but 0/i-events are enhanced. This 
ratio can be written as |25| 



~ jdN^y (10) 

The number of expected 0/t-events can be calculated as 

^ Ev) = [ ^k m + lt {K)P ^ ~* Ut) n{Ev) (1 " Br(T ^ (U) 

where we again supposed no contamination and the trigger and identification ef- 
ficiencies to be 100%. We can infer dN nc /dE u , the expected energy spectrum for 
nc-events in the MINOS far detector in the case of no flavor change, using the 
approximation 
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dN n 
dE L 



dN cc a 



cc \E V 



(12) 



The cross sections a h 



<j v -cc were taken from Ref. 26] and a 



from Ref. 



m 



In Table [TI|, we show our estimation of this ratio for no v ^ — > v T conversion, 
mass induced oscillation for several values of Am 2 and PDM with 7 = j^ est for 
the medium and high energy beam. For the low energy beam this test is ineffective 
since, in this case, MINOS will work essentially like a disappearance experiment. 
For the medium and high energy beams, the ratio predicted for the best-fit values 
of the parameters of the two solutions seem to be well separated. Nevertheless one 
has to take these results carefully, in a real experimental situation, experimental 
efficiencies, event contamination, resolution and systematic errors can substancially 
affect 0/x-event observables. 



C. CERN-to-Gran Sasso Neutrino Oscillation Experiments 

The new CERN neutrino beam to Gran Sasso is a facility that will direct a 
beam to the Gran Sasso Laboratory in Italy, at 732 km from CERN. Such a 
beam together with the massive detectors ICANOE/OPERA |22| in Gran Sasso 
will constitute a powerful tool for long baseline neutrino oscillation searches. The 
number of protons on target is expected to be 4.5 xlO 19 per year, the beam will 
have an average energy of 17 GeV, the fractions v e j v^fv^ and v T j in the beam 
are expected to be as low as 0.8%, 2% and 10~ 7 , respectively P7| . The number of 
charged current events per protons on target and kton, without neutrino oscillation, 
is calculated to be 4.7 x 10~ 17 , for 1 GeV < E v < 30 GeV |27|. These experiments 



are supposed to start taking data around 2005 and may be used to try to distinguish 
the two ANP solutions. 

We have investigated the capability of these detectors, in particular OPERA, 
working in the — > v T appearance mode, to elucidate which is the correct solution 
to the ANP. We have obtained the expected number of v T events, N T , that will be 
measured by OPERA, considering a pure beam, using the following expression 

N t = aJ (j)^{E) x P( Vft -> v T ) x Mr i) x <t„ t -cc(£) X e^E) dE, (13) 

i=l,h 

where <p u is the flux of at the Gran Sasso detector and <j Vt -cc is the charged 



current cross section for u T , both taken from Ref. [26]. The number of active targets 
A can be calculated as A = C n x M d x N A x 10 9 x N p x N y where M d is the detector 
mass in kton, x 10 9 is the number of nucleons per kton (where = 6.02 x 10 23 
is the Avogadro's number), N y is the number of years of data taking and N p is the 
number of protons on target per year and C n = 0.879 is normalization constant that 
we need to introduce in order to be able to reproduce the numbers presented in 
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Table 27 of Ref. []13||. We have assumed that v T identification will be accomplished 



through its one-prong decays into leptons (/) and hadrons (h). We have used the 
overall value J2i=i,h Br(r — > i) x £i(E) = 8.7 %, admitting a total mass of 2 kton for 



5 years of exposure, according to Ref. |T3[. The probability Piy^ — > v T ) is either 



supposed to be equal to zero (for no flavor transformation), to the MIO probability 
(Eq. (D) or to the PDM probability (Eq. (§)). 

In Table |T|, we show the number of u T events we have calculated, according 
to Eq. (0), for 5 year exposure, assuming sin 2 29 = 1 and four different values of 
Am 2 and the best-fit point of the decoherence solution to the ANP. We also quote 
in this table the total number of background events, which remains after all the 
kinematical cuts have been applied, normalized to an exposure of 5 years, taken 
from Ref. |13| . We observe that for the decoherence effect the rate of tau events is 
substantially higher than that of most of the Am 2 hypotheses and than the number 
of background events. 

We show in Fig. |7| the calculated energy spectrum for both flavor changing sce- 
narios at OPERA for 5 years of exposure, this was calculated using Eq. (|T3|). In 
spite of the fact that in some cases we observe that the ratio between the number 
of v T events coming from PDM and MIO is higher than four, in practice, it seems 
to be very difficult to observe this difference due to low statistics. 

In Fig. |8|, we show the result of the first statistical significance test. From this 
we see that discrimination between PDM and MIO can become difficult if Am 2 > 
4.4 x 10" 3 eV 2 {n a < 5). 

It is worth mentioning that if we vary the decoherence parameter in the range 
0.25 x 10" 21 GeV < 70 < 1.1 x 10~ 21 GeV, the number of expected N T events also 
varies from 31 to 114. This suggests that there might be some situations, even if 
the measured number of events is consistent with Am 2 < 4.4 x 10~ 3 eV 2 , where it 
might be difficult to disentangle the two ANP solutions in OPERA. This becomes 
clearer in Fig. ^, where we see that for data compatible with Am 2 < 4.1 x 10 -3 eV 2 
a separation of more than 3 a can be achieved for 70 > 4 x 10~ 21 GeV, but data 
consistent with values of Am 2 ~ 5.3 x 10~ 3 eV 2 can only be separated from the 
PDM solution if 7o > 7 x 10~ 21 GeV (n a > 5). 



D. Neutrino Factory in Muon Collider 



Many authors [ |T4| , pB| , pU | have emphasized the advantages of using the straight 
section of a high intensity muon storage ring to make a neutrino factory. The muons 
(anti-muons) accelerated to an energy (E^) constitute a pure source of both 
ip^) and v e {y e ) through their decay — > e~u e v^ (fi + e + z/ e z/ M ) with well known 
initial flux and energy distribution. 

In the many propositions for this type of long baseline neutrino factory one can 
find in the literature, the stored muon (anti-muon) energy E^ (Ep) ranges from 10 
GeV to 250 GeV and the neutrino beam is directed towards a faraway detector at a 
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distance corresponding to an oscillation baseline L varying from 730 km to 10 000 
km. 

Here we have explored the neutrino factory as a disappearance — > exper- 
iment. We will explicitly discuss the case when negative muons are stored in the 
ring, a similar calculation can be performed when positive muons are the ones which 
decay producing neutrinos. The relevant observable is the total number of /i that 
can be detected when \i is the produced charged lepton in the beam, i.e. the number 
of "same sign muons" that we will denote here by N^. We define N^, for unpolarized 
muons (see Appendix), as 



N, = n„M d — ^-f / h Q {x)e,{xE,) "_ K ^ [1 - -> u T )] 

+ G ^ XE ^ Br(r -> „) P(u, -> ds, (14) 

where is the muon source energy, x = E u /E^, M d is the detector mass in ktons, 
number of useful \i decays, 10 9 A^a is the number of nucleons in a kton and m M is 
the mass of the muon. The function ho contain the energy spectrum normalized to 
1 explicitly given in the Appendix. The charged current interaction cross sections 
per nucleon <j v _ cc and 0V T -cc can be found in Ref. PEfl . This number has two 
contributions, one from the produced in the decay that survive and arrive at the 
detector interacting with it producing a final /i, another from the v T that is produced 
by the flavor conversion mechanism and interact in the detector via charge current 
producing r, which subsequently decays to // with a branching ratio Br(r — > //). 
We have calculated the above integral from E th = 3 GeV to ensure good detection 
efficiency [^TJ, so that we can consider the \i detection efficiency to be 100 %, 



independent of energy. 

The muon beam is expected to have an average angular divergence of the 
O(0.1 m^/Efj). It was pointed out in Ref. that this effect is about 10 % so 
we have multiplied our calculated number of events by 0.9 to account for this. The 
contribution of the background to the number of muons observed in the detector, 
that includes muons from charm decays produced by charged current and neutral 
current interactions in the detector, has been neglected on account that this would 
be a small global contribution. 

Since the final configuration of such a facility is still not defined, we have tried 
to estimate the optimum configuration in order to maximize the difference between 
mass induced — > survival probability and pure decoherence. As observed in 
Ref. |J, at the best-fit point of the MIO and PDM solutions, the argument of the 
cosine, in the MIO case, and the argument of the exponential, in the PDM case, can 
be viewed to be approximately the same. Explicitly, 

P(u, -> Vp ) ~ + cos(A)] (MIO) , (15) 
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p(^^)~l[l + e -*] (PDM) , (16) 

where 

X = /3^, (17) 

with (3 = 2.54 Am 2 ~ 1.0 x 10 19 7 , Am 2 given in eV 2 , 70 in GeV and f3 in GeV/km. 
So that we can maximize the difference between the probabilities simply by finding 
the values of X, which maximize the function 

F(X) = \cos(X) -e~ x \. (18) 

This means that one has to numerically solve for X the transcendental equation 

e- x - sin(X) = 0, (19) 

with the condition cos(X) > — e~ x . 

Once we find the spectrum of solutions X we can apply Eq. ( |T7D to find the 
optimal distance L opt by choosing a value for and fixing E v at the average value 
of the observable neutrino energy, i.e.. (E u ) = 0.7 E^ |3(|, with /3 = 6.1 X 10~ 3 



GeV/km, the best-fit value of Ref. ||. We have calculated the maximal difference 
between the two survival probabilities for E^ = 10 GeV, 20 GeV, 30 GeV and 50 
GeV as a function of L, but only kept the cases where this difference reaches about 
50 %. This was obtained for L opt = 3096 km, 6192 km, 9289 km and E^ = 10 GeV, 
20 GeV and 30 GeV, respectively. 

Note that using this criteria, we are looking for a maximum in the absolute dif- 
ference between the survival probabilities ((9(0.5)), which is a very strong condition 
that we impose in order to detect and distinguish the signal of disappearance. 
Since this depends on the number of muon decays, which is big, this ensures that 
the signal will be quite sizable experimentally. We have seen before (see for instance 
Fig.0), that there are cases of L/E where the absolute difference in the survival 
probabilities is small while the ratio Piy^ — ► u T )(PDM)/P(u tl — ► u T )(MIO) is big, 
but since the statistics is poor this ratio is an ambiguous observable to establish the 
difference between PDM and MIO. 

Using Eq. (|i"4]) we have estimated the neutrino event rate, for the three different 
hypotheses: no flavor conversion, MIO and PDM. 

We addressed a few optimal configurations so that our conclusions may be useful 
in planning real experimental setups. For each optimal situation, we have performed 
our calculation for two possible scenarios: (1) a total of 1.6 x 10 20 /1 decays per year 
with a detector of 10 kton; (2) a total of 2.0 x 10 20 \i decays per year with a detector 
of 40 kton. Since the event rate is proportional to E^/L 2 (see Eq. flT4"D), the most 
conservative configuration is the one with L opt = 3096 km and E^ = 10 GeV. Our 
results are summarized in Table [IV]. It is clear that one can distinguish among all 
the studied hypotheses, even in the most modest scenario (1). 
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An even more conclusive result can be obtained by looking at the energy dis- 
tribution dN^/dE u (see Eq. ( |A10| )) for charged current \i events. We show, in Fig. 
10 , our predictions, for five years of data taking, assuming 2.0 xlO 20 /i decays per 



year, = 10 GeV, L = 3096 km for null oscillation, the best-fit point of the PDM 
solution to the ANP and four different values of Am 2 which are consistent with the 
MIO solution to the ANP. We see that in the majority of distinct signal will 

be observed, making it possible to establish which of the flavor changing hypotheses 
is the correct one. 

Although we have optimized the setup parameters in order to maximize the 
separation between the two ANP solutions at their best-fit point, we also have 
checked that even if we allow 70 to vary around the best-fit value, i.e. 0.25 x 
10 -21 GeV < 70 < 1.1 x 10~ 21 GeV, we will get an energy spectrum which will 
always increase monotonically with the neutrino energy, in this range these curves 
will all coinciding up to 4 GeV, after that if 70 < 0.6 x 10~ 21 GeV the growth will 
be slightly slower than what is shown in Fig. [H], for 70 = 0.6 x 10~ 21 GeV and if 
70 > 0.6 x 10~ 21 GeV the curves will be steeper approaching the one for no flavor 
change. 

In Fig. 11, we show the statistical significance test for 5 years of data taking in the 
plane Am 2 x sin 2 29 assuming the most conservative of the optimal configurations. 
From this figure, we see that for most of the MIO parameter space the two solutions 
can be well discriminated. This conclusion is further confirmed by Fig. 12, where 
the statistical significance test was performed letting 70 vary. Here we note that 
for data compatible with 1.5 x 10~ 3 eV 2 < Am 2 < 4.8 x 10 -3 eV 2 a very clear 
separation is possible V70. It seems that unless the data prefers lower values of 
Am 2 < 1.5 x 10~ 3 eV 2 or a few isolated islands in the plane 70 x Am 2 , in which case 
a clear separation between the two solutions might be compromised, decoherence 
and mass induced oscillation will present a very distinct signature of their dynamics. 
Fortunately, for Am 2 < 1.5 x 10~ 3 eV 2 and for Am 2 > 5.0 x 10~ 3 eV 2 a clear cut 
between the two solutions can be accomplished respectively by the OPERA and the 
MINOS (medium) experiments, as we have discussed in Sec. [IV C| and [IV B. 



V. CONCLUSIONS 

We have discussed the perspectives of future experiments of distinguishing the 
MIO solution to the ANP from the PDM one. This is specially important since 
it will be very difficult for NOMAD and CHORUS to achieve a sensitivity in the 
Vfj, — > v T mode to directly exclude/confirm the latter solution [17] . 



Our study on discriminating these two ANP solutions permitted us to arrive at 
the following general conclusions: K2K and probably MINOS will not be able to 
shed much light on the dynamics which promotes — > v T conversion, OPERA and 
a neutrino factory in a muon collider are more suitable for this job. 

From the statistical significance tests we have performed, considering positive 
discrimination only if n a > 5, we can say the following. K2K cannot discriminate 
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PDM from MIO if the data is compatible with 2.2 x l(T 3 eV 2 < Am 2 < 4.5 x l(T 3 eV 2 , 
V70. In fact, this almost covers the entire parameter space allowed by the SuperK 
atmospheric data at 99% C.L., as can be seen in Fig. |2|. MINOS, in the low energy 
beam configuration, cannot separate the two solutions if the data is compatible 
with 1.1 x 10~ 3 eV 2 < Am 2 < 5.2 x 10~ 3 eV 2 , V70; in the medium energy beam 
configuration if 3.0 x lCT 3 eV 2 < Am 2 < 4.0 x 10~ 3 eV 2 , V70; in the high energy 
beam configuration if 4.8 x 10~ 3 eV 2 < Am 2 < 5.1 x 10~ 3 eV 2 , V70. Running with 
the low energy setup, MINOS will be very similar to K2K, and certainly will not 
be able do discriminate the solutions, see Figs. [5] and |6|. For the high energy setup 
MINOS will be more selective and similar to OPERA. OPERA, after 5 years, will 
only not disentangle PDM from MIO if the data is compatible with 5.9 x 10~ 3 eV 2 < 
Am 2 < 6.7 x 10~ 3 eV 2 , V70. This corresponds to the upper corner of the SuperK 99% 
OL. region, see Fig. [8]. A neutrino factory with = 10 GeV, L = 3096 km and 1.6 
xlO 20 fi~ decays per year, with a detector of 10 kton after 5 years of data taking will 
still not be able to discriminate the two solutions if the data prefers a small island in 
the plane sin 2 29 x Am 2 , at Am 2 ~ 1.5 x 10~ 3 eV 2 and sin 2 29 ~ 0.95, see Fig. 0. On 



the other hand, if the data is compatible with 1.5 x 10~ 3 eV 2 < Am 2 < 5 x 10 _3 eV 2 
the separation between PDM and MIO will be extremely clear, V70. This type 
of facility, it seems, will be the only one able to measure spectral distortions and, 
thereby, directly test if neutrino flavor change is indeed an oscillation phenomenon. 

In any case, the combination of the results of all these proposed experiments will 
most certainly unravel the dynamics of neutrino conversion in — > v T mode, if only 
active neutrinos exist in nature. 

There is a proposed atmospheric neutrino experiment, MONOLITH ||32| , which, 
in principle, could distinguish PDM from MIO through the observation of the first 
oscillation minimum. Here we have only investigated the capabilities of accelerator 
neutrino experiments which we believe have to be performed in order to completely 
ratify the mechanism behind neutrino flavor change. 
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TABLE I. Expected number of i^-events calculated for the K2K Far Detector of 22.5 
kton. In all cases sin 2 29 = I. 



Hypothesis 


iV FD 


Number of POT 


No Flavor Change 


40.3 


2.29 x 10 19 


Am 2 = 3 x 1(T 3 eV 2 


28.3 




Am 2 = 5 x 1(T 3 eV 2 


18.9 




Am 2 = 7 x 1(T 3 eV 2 


13.3 




Am 2 = 3 x 1(T 3 eV 2 


123 


1.0 x 10 20 


7o = 0.6 x 10" 21 GeV 


122 
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TABLE II. Expected ratio Rq^/i^ in MINOS for 10 kton year exposure. The error is 
only statistical. In all cases sin 2 26 = 1. 



Medium Energy Beam 



Hypothesis 




No Flavor Change 




0.33 ±0.01 


Am 2 = 1.5 x 1(T 3 eV 2 




0.35 ±0.01 


Am 2 = 3.0 x 10~ 3 eV 2 




0.41 ±0.01 


Am 2 = 4.5 x 10~ 3 eV 2 




0.52 ±0.01 


Am 2 = 6.0 x 10 -3 eV 2 




0.71 ±0.01 


70 = 0.6 x 10" 21 GeV 


0.47 ±0.01 


High Energy Beam 


Hypothesis 


Ron/in 


No Flavor Change 




0.314 ± 0.004 


Am 2 = 1.5 x 10" 3 eV 2 




0.320 ± 0.004 


Am 2 = 3.0 x 10~ 3 eV 2 




0.340 ± 0.004 


Am 2 = 4.5 x 10" 3 eV 2 




0.380 ± 0.005 


Am 2 = 6.0 x 10" 3 eV 2 




0.430 ± 0.005 


70 = 0.6 x 10" 21 GeV 


0.413 ±0.005 



TABLE III. Expected number of v T -events calculated for OPERA assuming 5 year 
exposure, using 8.7% for the total r selection efficiency. In all MIO cases sin 2 20 = 1. 



Hypothesis 




Am 2 


= 1.5 x 10~ 3 eV 2 




4 


Am 2 


= 3.0 x 10~ 3 eV 2 




16 


Am 2 


= 4.5 x 10~ 3 eV 2 




36 


Am 2 


= 6.0 x 10" 3 eV 2 




62 


7o = 


0.6 x 10~ 21 GeV 


69 


Total Background ~ 0.6 events 
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TABLE IV. Expected number of fi events, N^, for a few configurations of a neutrino 
factory under the three studied hypotheses. Scenario (1) corresponds to a total of 1.6 x 10 20 
H decays per year and a 10 kton detector, scenario (2) to a total of 2.0 x 10 20 fi decays 
per year and a 40 kton detector. The numbers in scenario (1) are calculated for one year 
while in scenario (2) are for five years of data taking. The errors given are only statistical. 
For MIO we used sin 2 26 = 1. 



Hypothesis 


(1) 


(2) 


L (km) 
Ey. (GeV) 


3096 
10 


6192 
20 


9289 
30 


3096 
10 


6192 
20 


9289 
30 


No Flavor Change 


1315 ± 36 


2580 ± 50 


3804 ± 61 


32880 ± 181 


64501 ± 253 


95101 ± 308 


Am 2 = 3 x 10" 3 eV 2 


222 ± 15 


503 ± 22 


805 ± 28 


5554 ± 75 


12596 ± 112 


20141 ± 141 


70 = 0.6 x 10- 21 GeV 


716 ± 27 


1439 ± 37 


2156 ±46 


17918 ± 134 


35991 ± 189 


53910 ± 232 
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FIG. 1. Spectral distortion expected for one ring /x-like events at K2K assuming 
— > f T flavor conversion with sin 2 2f9 = 1 and different values of Am 2 (MIO) and 

70 = 0.6 x 10~ 21 GeV (PDM). The bars represent the statistical error at each bin. The 

calculation was done for a total of 10 20 POT. 
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sin 2 20 

FIG. 2. Regions in the sin 2 26 x Am 2 plane, for 70 = 7o est , where the number of ./Vfd 
events expected for PDM and MIO are separated by u a = 1,2,3 and 5 for K2K after 
10 20 POT. The inner part of the gray curve is the one allowed at 99 % C.L. by the latest 
SuperK atmospheric neutrino data || . The sensitivity of K2K is marked by an horizontal 
line with an arrow. 
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Am 2 (eV 2 ) 



FIG. 3. Regions in the 70 x Am 2 plane, with sin 2 29 = 1, where the number of A^fd 
events expected for PDM and MIO are separated by n a = 1,2,3 and 5 for K2K after 10 20 
POT. The dotted gray lines mark the region allowed at 99 % C.L. by the latest SuperK 
atmospheric neutrino data || and the cross the best fit values of the PDM (8| and the 
MIO H solutions to the ANP. The start of the sensitivity of K2K is marked by an arrow. 
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FIG. 4. Spectral distortion expected at MINOS for the best-fit points of the vacuum 
oscillation and of the pure decoherence solutions to the ANP for three possible beam 
configurations. 
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FIG. 5. Same as Fig. ^ but for Ni„ events in MINOS after 10 kton-year exposure. 
The sensitivity of MINOS is marked by an horizontal line with an arrow. 
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FIG. 6. Same as Fig. || but for iVia events in MINOS after 10 kton-year exposure. The 
sensitivity of MINOS is marked by a vertical continuous line with an arrow. 
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FIG. 7. Spectral distortion expected at OPERA for different values of Am 2 (dashed 
lines) as well as for the best-fit point of the pure decoherence (continuous line) solution to 
the ANP, for 5 year exposure. We have taken into account the overall efficiency of 8.7 % 
in accordance to Ref. |Ol. In all cases sin 2 20 = 1. 
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FIG. 9. Same as Fig. || but for iV T events in OPERA after 5 years. The sensitivity of 
OPERA is marked by an arrow. 
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2.0 X 10 jll decays per year 




E„(GeV) 

FIG. 10. Spectral distortion expected at a neutrino factory for = 10 GeV, L = 3096 
km and several values of Am 2 as well as for no flavor change (continuous line) and for 
the best-fit point of the pure decoherence (dot-dashed line) solution to the ANP, after five 
years of data taking. In all MIO cases sin 2 29 = 1. 
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FIG. 11. Same as Fig. |2| but for jV„ events after 5 years of a neutrino factory in the 
scenario (1) of Table IV. The sensitivity of this setup is marked by an horizontal line with 
an arrow. 
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FIG. 12. Same as Fig. ||| but for Nu events after 5 years of a neutrino factory in the 



scenario (1) of Table IV. The start of the sensitivity of this setup is marked by an arrow. 
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APPENDIX: 



The distribution of in the decay fi — > e + v e + in the muon rest-frame 



Jem) is given by |L4| 

d 2 N Vii 1 



( j r(J Q —[h (x)+V ll h 1 (x)cose cm } } (Al) 

x = 2E£ m /m /J , where denotes the neutrino energy, 8 cm is the angle between 
the neutrino momentum vector and the muon spin direction, and is the average 
muon polarization along the beam directions. The functions h and hi are given in 
Table 0. 

On applying a Lorentz transformation to boost into the laboratory frame (lab), 
it is found that the neutrino energy distribution at any polar angle is just scaled by 
a relativistic boost factor depending on the angle 



1 

2 " 

and 



^ ab = -x4 ab (l + /?cosM , (A2) 



• a sin# cm 

sin #i ab = . 7j T > A3 

7(l + pcos# c ' 



' cm/ 



where (3 and 7 are the usual relativistic factors and we have used 7 = Erf jm^. 

Because in the long-baseline experiments only emitted in the forward direc- 
tion are relevant to the computed flux, one can make the following approximation: 
cos^cm ~ 1, sin# cm ~ 9 cm . Also at high energy /3 ~ 1. Hence Eq. (|A2[) leads to 
x = E^/E 1 ^ . Then, we can rewrite df2 cm in terms of df2i a b as 

rffi cm = 7 2 (l + /5) 2 dn lab . (A4) 
Substituting Eq. (|A4|) into Eq. (|AT|), we obtain as a function of the lab variables 

-T-^- = ( x ) ± P m ^1 (*) cos £ cm ] , (A5) 

which for unpolarized muons simplify to 

d 2 N^ {E^f dN u ^ E^ 

= ho {x) -^dE^ = 47r ^ ho{x) ■ (A6) 

The number of expected /1 events from interactions with the detector is given 

by 

JV /1 En (1 M d 10 9 iV A x$x ( 7Xe fl (A7) 

V ■* 

normalization 
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The flux $ can be written as 



N v .. 1 



1 dN v 



AttL 2 t dE l f° AttL 2 t dE l * h ' 

Substituting Eq. ( POp into Eq. (|A8|), we get 



(A8) 



-h Q (x). 



m 2 irL 2 t 



(A9) 
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TABLE V. Flux functions h (x) and h^x). 



h (x) 


hi(x) 


2x 2 (3-2x) 


2x 2 (l-2x) 



Finally, after t-years 

dN, _n,M d 10 9 N A _ lab 



E l ™h (x) x a x e M . (A10) 
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